090  id’s 


DEGENERACY  IN  IDEAL  CHEMICAL 
EQUILIBRIUM  PROBLEMS 

BY 

JAMES  H.  BIGELOW 

TECHNICAL  REPORT  NO.  70-3 
MARCH  1970 


D  C 

W  MAY  21  :2--o  ij 

jurasaTiu 

E 


TfU&  c'xfcvment  hat  been  appfwv^d  pubtCc.  nztvut  and  iote.; 
aJU  yUAtfubution  lA  antmitid. 


Operations 

Research 


House 


\ 


Stanford 

University 


CALI  FOR  NI A 


DEGENERACY  IN  IDEAL  CHEMICAL 
EQUILIBRIUM  PROBLEMS 


by 

JAMES  H.  BIGELOW 


TECHNICAL  REPORT  NO.  70-3 
MARCH  1970 


Department  of  Operations  Research 
Stanford  University 
Stanford,  California 


Research  and  reproduction  of  this  report  was  partially  supported  by 
Office  of  Naval  Research,  Contract  N-00014-'67-A-01 12-0011;  U.S.  Atomic 
Energy  Commission,  Contract  AT(04-3}  326  FA#18;  National  Science 
Foundation,  Grant  GP  9329;  U.S.  Army  Research  Office,  Contract 
DAHC04-67-nn?8;  and  National  Institute  of  Health,  Grant  GM  14789-02. 


Reproduction  in  whole  or  in  part  is  permitted  for  any  purpose  of  the 
United  States  Government.  This  document  has  been  approved  for  public 
release  and  sale;  its  distribution  is  unlimited. 


CONTENTS 


Section 

I.  INTRODUCTION . 1 

II.  VIRTUAL  MOLE  FRACTIONS . 5 

III.  OVERCOMING  DEGENERACY  WITH  SLACKS . 20 

1.  Slacks:  General  Form . 20 

2.  Slacks:  Special  Form . 27 

REFERENCES . 48 


DEGENERACY  IN  IDEAL  CHEMICAL 
EQUILIBRIUM  PROBLEMS 
James  K.  Bigelow^ 

The  RAND  Corporation^  Santa  Monica,  California 


I,  INTRODUCTION 

>  A  single-  or  multi-phase  chemical  equilibrium  prpbl^ 
may  be  thought  of  as  the  problem  of  minimizing  a  particular 
nonlinear  function  (the  free  energy)  of  composition  subject 
to  the  conditions  that  the  composition  vector  be  nonnegative 
and  satisfy  a  system  of  linear  equations 
laws) .  It  was  pointed  out  in  a  previous  ; 
free  energy  is  convex  and  homogeneous  of  degree  one,  but 
that  as  a  variable  approaches  zero,  the  free  energy  may 
behave  badly. 

^his  P£pe5,  the  second  in  a  series  of  three**, 

—/the  phrase  "chemical  equilibrium  problem""^refers  only  to 
a  problem  V7ith  a  particular  mathematical  form.  Problems 
of  this  form  arise  in  many  situations  that  are  not  clas¬ 
sically  denoted  chemical  equilibrium  problems For  example, 

*Any  views  expressed  in  this  paper  are  those  oX  the 
author.^  They  should  not  be  interpreted  as  reflecting  the 
views  of^The  R(\K3D  Corporation  or  the  official  opinion  or 
policy  of  any  of  its  governmental  or  private  research  sponsors, 
Papc3.s  are  reproduced  by  Tne  Corporation  as  a  courtesy 

to  members  of  its  staff. 


(i:h«Hnauu43filAnce 


**The  first  is  ‘'Chemistry,  Kinetics  and 
the  third,  ‘ Computing  Equilibrium  Composition 
Chemical  Systems". 
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the  dual  to  a  geometric  programming  problem  [2],  [3],  has 
this  form.  Also^  steady— state  problems  [4]#  ma'iy  of  which 
arise  naturally  in  industry  and  in  the  chemical  laboratory, 
can  often  be  represented  mathematically  by  problems  of 
this  form. 

We  call  the  problems  under  (Consideration  "ideal”  be¬ 
cause  they  are  derived  from  chemical  equilibrium  problems 
which  take  the  simplest  form  possible.  As  a  consequence, 
the  free-energy  function  is  as  simple  as  it  can  be  [1], 

The  ideal  chemical  equilibrium  problem,  then,  is  the 
problem  of  minimizing  a  function  F(xj^,  X2,  j 

finded  below,  subject  to  the  linear  constraints: 

n 

(1.1)  S  a. .x.^b.  ,i— 1,  2,  ...,m 

and  the  inequalities: 

(1.2)  j  j  =  jL,  2,  ...,  n. 

The  a.,  and  b.  are  given  real  constants.  We  assume  that 
the  m  equations  (I.l)  are  linearly  independent,  and  so 
that  the  problem  will  not  be  trivial,  that  m  <  n. 

The  n  variables  Xj  arc  partitioned  into  p  noncn.pty 
subsets  called  coaipartmcnts,  or  phases.  We  denote  the 
compartment  containing  the  j — ^  variable  by  <j>.  Wo  may 
indicate  that  x^  and  X|i^  arc  in  the  same  compartment  by 
writing: 
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j  e  <k>  or  k  e  <j>  or  <j>  “  <k>. 


Sach  compartment  has  associated  with  It  a  sum. 


(1.3) 


S  X.. 
J6<k>  ^ 


Each  variable  has  associated  with  it  a  variable  fraction. 


(1.4) 


The  objective  function  to  be  minimized  Is! 


(1.5) 


^  ^^1^  •  •  •  # 


n 

E 

J»1 


+  log  i^j) . 


The  quantities  c^,  c^,  ...,  c^  are  given  real  constants. 

If  Xj  0  for  some  j,  then  either  Xj  *  0  (in  the  event 
that  >  0)  or  Xj  is  undefined  (if  0) .  In  either 

case,  log  x^  is  undefined;  but  to  maintain  the  continuity 
of  F  at  the  boundary  of  the  constraint  set,  we  define 
Xj  log  Xj  «  0  whenever  x^  «  0  ([5]>  p.  364). 

It  will  be  convenient  to  use  matrix  notation.  We 
let  A  be  a  matrix  whose  ij~  element  is  (The  first 

index  refers  to  the  row  number,  the  second  to  the  column. 
Thus  A  has  rn  rows  and  n  columns.)  Similarly,  b  is  a 
column  vector  of  dimension  ra,  whose  1 —  component  is  b^. 


In  the  same  fashionj  we  let  x  be  the  n~vector  with 
coinponents  x^ ;  x  the  n-vector  with  components  Xj ;  log  x 
the  n-vector  with  components  log  x^;  and  c  the  n~vector 
with  components  Cj,  All  of  these  are  column  vectors. 

Using  vector  notation,  we  may  wite  the  chemical 
equilibrivim  problem  in  the  follo\dLng  compact  form: 

Min  F(x)  =  Min  x  •  (c'  +  log  x) 

(1.6)  s.t.  A2c  =  b 

X  >  0. 

The  notation  x  •  (c  +  log  x)  denotes  the  inner  product 
of  t:he  vectors  x  and  (c  -h  log  >:) . 

As  we  shall  see,  the  singular  behavior  of  F  at  the 
boundary  of  its  domain,  the  ncmegative  orthant,  makes  it 
difficult  even  to  recognize  a  solution  to  (1.6).  In  this 
paper  we  discuss  a  method  for  surmounting  this  problem, 
called  degeneracy,  and  a  method  for  avoiding  it. 
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II.  VIRTUAL  MOLE  FRACTIONS 

In  this  section  we  will  discover  how  to  recognize 
when  a  vector  x  is  a  solution  to  ProhX^  (1.6).  We 
touched  on  this  question  in  [1]^  where  we  showed  a  way 
to  recognize  whether  x  solved  (1.6)  if  x  were  strictly 
positive.  In  this  section  we  extend  the  result  to  include 
all  feasible  x. 

It  will  be  convenient  to  define  the  concept  of  an 
admissible  direction. 

Definition;  Let  a  cranposition  vector  x  be  feasible  in  a 
given  chemical  equilibrium  problem.  That  is,  x  satisfies 
Ax  “  b,  X  >  0,  An  n-vector  6  is  said  to  be  admissible 
at  X  if  for  every  t  >  0,  t  sufficiently  small,  x  *1  te  is 
feasible. 

We  will  have  occasion  to  consider  several  different 
chemical  equilibrium  problems  with  different  constraint 
sets.  When  we  refer  to  the  concept  of  admissibility,  we 
mean  that  a  direction  0  is  admissible  for  the  problem  then 
under  consideration. 

Mote  that  0  is  admissible  at  x  if  and  only  if  A0  =>  0 
and  0^  >  0  for  every  j  such  that  x.  “  0, 

J  J 

The  difficulty  in  recognizing  solutions  arises  from 
the  behavior  of  the  Gibbs  function  F(x)  when  some  of  the 


X.  <=  0.  It  is  well  known  that  since  F(x)  is  convex,  F 
achieves  a  minimum  ac  x  if  and  only  if  for  every  reaction 
vector  0  admissible  at  x, 

Fj(x)  >  0 

where  F^x)  is  the  derivative  of  F  at  x  in  the  direction 
6.  This  derivative  is  linear  in  G  whenever  x  >  0;  but  it 
may  be  non-linear— indeed  it  may  be  infinite — when  some  x.  «  0, 

Much  of  the  content  of  this  chapter  was  anticipated  by 
N.  Z.  Shapiro,  who  proved  a  restricted  form  of  the 
main  theorem,  II. 4.  Theorem  II. 4  in  all  its  generality  was 
proved  by  the  late  Dr.  Jon  Folkman,  although  the  present 
author  never  sav7  that  proof.  The  proof  appearing  here  is 
original  with  the  author. 

It  will  be  convenient  at  this  point  to  define  two  sets. 

H(A,  b)  =  {x  e  e’^IAx  =  b,  >  0) 

M(f1h(A,  b))  =  {x  e  H(A,  b)|F(x)  <  F(y) V  y  e  H(A,  b)}. 

(One  might  as  easily  define  M(F|W)  for  a  general  set  and 
general  function  F:W  -  E^.  Thus  M(F|VJ)  =  [x  e  W|F(x)  <  F(y) 

V  y  e  W}.)  Then  H(A,  b)  is  the  set  of  all  x  which  satisfy  both 
the  mass  balance  and  non-negativity  constraints  of  problem 
(1.6),  i.e.,  H(A,b)  is  the  sot  of  feasible  compositions  of 
the  chemical  system.  M(f|h(A,  b))  (written  M(f|iI)  if  there 
will  be  no  confusion)  is  the  set  of  solutions  to  problem 
(1.6),  and  hence  the  set  of  equilibrium  compositions  of 
the  cliemical  system. 
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In  this  section  we  will  derive  a  characterization  of 
The  name,  ’’Virtual  Mole  Fractions,”  arises  from 
the  fact  that  in  equilibrium  chemistry  it  is  convenient 
to  express  the  variables  in  units  of  ’’moles”  and  to  call 
the  variable  fractions  ’’mole  fractions.”  It  may  happen 
that  at  equilibrium,  one  or  more  phases  of  a  chemical 
system  are  empty.  In  this  case,  the  mole  fractions  of 
the  variables  in  the  vanished  phases  are  undefined;  but 
we  will  show  how  to  define  quantities  to  take  their  place, 
quantities  called  ’’virtual  mole  fractions.” 

Before  V7e  can  do  any  of  the  real  work  of  this  section 
we  must  satisfy  ourselves  that  Problem  (1.6)  may  be  solved 
using  Lagrange  multipliers. 

Theorem  II. 1:  Assume  there  exists  x  e  H(A,  b)  such  that 
X  >  0,  and  suppose  that  M(f1h(A,  b))  is  nonempty.  Let 
R  =  {x  £  E’’'^jx  >  0}.  Then  there  exists  tt°  e  e’^  such  that: 

M(f1r(A,  b))  c  M(F(x)  -  n°(Ax  -  b)|R). 

Ro;:.ark :  This  statement  says  exactly  that  each  solution 
of  (1.6)  must  also  minimize  the  Lagrangian  form 
F(x)  -  "^(Ax  -  b)  cimong  all  non-negative  x. 

Proof ;  From  [5],  Theorem  8.13,  p,  368,  we  know  that  F  is 
convex  on  R.  Hence,  the  Lagrangian  is  convex. 

Let  x°  e  M(fIh(A,  b))  (which  we  assumed  was  nonempty), 
and  define  the  compact,  convex  set  ~  >  Oj  ;'|X°  —  <  1] . 
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From  [6],  Theorem  6,  p,  478,  we  know  that  there 
exists  e  such  that.’ 

(11.1)  F(x°)  *=  Min  [F(x)  -  tt®(Ax  -  b)], 

xc 

*1  1 

Suppose  FCx"^)  <  F(x°)  for  some  x"^  e  R.  By  choosing  1  >  X  >  0 
sufficiently  small  we  may  construct  y  =  (Xx^  +  (1  —  X)x°)  e 
which  by  convexity  must  satisfy  F(y)  <  F(x°).  Thus: 

(11.2)  F(x°)  =  Min  [F(x)  -  Tr°(Ax  -  b)]. 

xeR 

The  result  is  immediate.  Q.E.D. 

For  any  x  >  0  we  partition  the  indices  fl,  2,  n} 

into  two  disjoint  sots: 

1  {j \^<j>  “ 

d  =  Cj  >  Oj. 

Given  any  vector  9  g  E^j  we  separate  it  into  two  terms, 
and  6j, 

(  G.  if  j  e  I 

(®T)i  = 

1 0  othen>7ise 


With  the  above  notation,  we  may  use  [5],  Theorems 
8.10  and  8.11,  p.  368,  and  the  vector  notation  developed 
in  section  I  to  write  the  derivative  of  F  at  x  in  the 
direction  8  in  the  followitig  compact  form; 
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(II.3)  F*(x)  =  ej»(cj  log  Sj)  +  0j*(cj  +  log  Xj), 

This  is  defined  for  all  0  admissible  at  x.  We  understand 

that  F*(x)  may  take  on  the  value  Thus  if  for  some 

j  €  J  we  have  x.  =  0,  we  see  that  log  x.  =  In  this 

J  ,  J 

case,  0.  >  0  implies  that  Fg(x)  ~  -®. 

•3  • 

Theorem  II.  2;  Let  x  6  H(Aj  b)  .  Then  x  e  M(F|H)  if  and 

I 

only  if  F  (x)  >  0  for  all  0  admissible  at  x. 

b 

Proof:  [5],  Theorem  10.2,  p.  372. 

QED. 

Corollary  II. 3:  Let  x  e  H(A,  b)  and  X  >  0.  Then  x  €  M(FIH) 
if  avid  only  if: 

(XI. 4)  8  •  (c  -r  log  x)  =  0  V  G  3  A0  =  0. 

Proof:  [5]j  Theorem  10.3,  p.  372.  This  also  follov7s 

imniodiately  from  Theorem  11,2,  since  x  >  0  implies  I.  is  i 

eioty. 

.  QED. 

Notice  that  if  x  e  H(A,  b),  and  x  >  0,  the  optimality 
condicion  in  Corollary  11,3  requires  a  kno\:lcdge  only  of  | 

I 

not  of  X.  It  is  the  mole  fractions,  not  the  moles  that  count.  1 

[ 

In  the  event  that  I  is  not  empty,  wo  :uld  like  something  f 

L 

to  take  the  place  of  x  in  (II. 4) .  This  motivates  the 
folloxv’ing  definition  and  theorem. 


9 


(II.  5a)  §(x)  «  {cp  c  E’^lcp->  0,  8*(c  +  log  cp)  =.0  ‘  3  As  =  0, 

and  cp.  =  v  j  c  J. 

3  3 

Eltiiientary  linear  algebra  shows  that  we  may  write  the 
equivalent  definition: 

(II.  5b)  §(x)  =  [cpeE’^Icp>0,  c  +  log  <p  =  A^tt  for  some  rr  e  E™, 

► 

and  <P;r  “  V  j  e  J. 

Theorem  II. 4:  Assume  S  y  e  H(Aj  b)  3  y  >  0,  and  let 

X  e  H(A,  b).  Then  x  e  M(f1h)  if  and  only  if 

a  cp  s  5(x)  3  <-l  for  every  <j>,  where  ~  «Pk. 

3  ■“  3  ke<j> 

To  prove  Theorem  II. 4 j  we  V7ill  need  the  follov/ing 
lemma. 


Lemma  II. 5:  Let  cp  e  satisfy  cd  >  0.  Then  for  any 
8  e  eP,  8^0, 


^  P  _  P  0^ 

Proof;  Let  8  =  L  81,^  P  =  L  m,  j  8 .  =  —I,  cb.  =  --3. 

k=l  k=l  ^^3  0  3  ^ 

Rearranging  terms  in  (II. 6)  we  find  it  equivalent  to: 

(11.7)  I  0.  (log  e.  -  log  cp.\  >  0. 

j=l  3  \  3  3/  - 

This  function  of  0  is  of  exactly  the  same  form  as  Gibbs 
function,  v?ith  c,  ••=  -log  rp.  If  we  minimise  the  function 

J  J 

in  (II. 7)  subject  only  to  the  constraint  0  >  0,  we  V7ill  be 
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solving  a  chemical  equilibritsn  problen.  The  absence  of 
mass  balance  constraints  tells  us  that  every  vector  v  e 
is  a  reaction  vector. 

Let-  6  >  0,  and  take  the  derivative  in  the  direction 
V.  By  <11. 3)  this  becomes: 

F^(0)  =>  E  VjClog  e.j  -  log 

Clearly  this  is  non-negative  for  every  v  e  E^  if  and  only 
if  for  each  j, 

0.  «=  cp.  , 

Thus  the  function  is  minimized  if  and  only  if  for  some 
a  >  0,  0  e  acp. 

But  clearly  e  -  implies  that  (II. 7)  is  zero.  Q.E.D. 

Proof  of  Theorem  II. 4:  Define  the  Lagrangian  problem  P(t7) 
to  be; 

(II. 8)  P(t7)  :  Min  G(x)  =  Min  [F(x)  —  n(Ax  —  b)  ] 

s.t,  X  ^  0. 

Let  x°  6  H(A,b)  .  By  Theorem  II. 1^  it  will  be  sufficient 
to  silo's?  that  the  following  statements  are  equivalent: 

(i)  S  cp  e  $(x°)  3  ^  1  V<j> 

(ii)  a  TT°  e  e’"'  3  x*^  solves  P(tt°). 
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Note  that  for  each  tt,  P(tt)  is  just  a  chemical  equi¬ 
librium  problem  with  no  linear  equality  constraints.  Thus 
by  Theorem  II. 2,  solves  P(n°)  if  and  only  if 

(II.  9)  Gq(x°)  =  +  log  Gj)  +  ^j)  ~  6*A^Tr°  >  0 

V0  admissible  at  x°. 

Clearly,  0  is  admissible  at  x°  in  this  problem  if 

and  only  if  0.  >  0  whenever  x?  =  0.  Immediately  V7e  see 
J  ~  j 

that  if  X?  ~  0  for  some  j  e  J,  v/e  can  find  8  with  Gi  >0, 

J  j  j 

and  force  G. (x°)  =  — Thus  we  may  suppose  x.  >  0  for  each 

W  J 

j.  e  J. 

It  is  obvious  that  (II. 9)  cannot  hold  uiilsss, 

(Il.lOj  Cj  log  xj  =  aJit®  . 

If  we  define  a  vector  c?  by: 

(IX.  11)  c  log  cp  A^7T°. 

then  clearly,  since  x.  >  0  for  j  e  J,  cp  e  5(x  ).  (See 

•J 

definition  (II. 5b).)  Substituting  (XI. 11)  into  the  inequality 
(II. 9),  we  find  that  x'^  solves  P(Tr°)  if  and  only  if: 

(11.12)  e-|-*  (log  -  log  cp^)  >0  Gj  >  0* 

Let  !■>:•  be  the  set  of  compartmciits  <j>  xdiose  indices 
are  in  I,  Then  (11.12)  becomes: 

(11.13)  D  (  Z)  e.(iog  0.  -  log  cp ))  >0  V  e_  >  0. 

<k>er>:  'je<k>  ^  ^  ^  ^  “ 
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Lemma  II. 5  considered  exactly  such  quantities  as  the 
inner  suun.  Thus  (11.13)  holds  if  and  only  if: 


-  D 

<k>e I* 


s' 


9j  1°8  'P<k>  S:  0 


V 


It  is  easy  to  see  that  this  is  true  V  e-j-  ^  0  if  and  only 
if 

(11.14)  ^1  V  <k>  e  I*  ■. 

Since  cp,  «  x  and  D  x.°  *=  1  if  k  g  J. 

jG<k>  J 

(11.14)  holds  if  and  only  if  statement  (i)  above  is  true, 
completing  the  proof.  Q.E.D. 


For  each  x  e  M(F|H),  we  ma3r  define  a  subset  of  §(x) 
as  follows: 


(11.15) 

We  call 
with  X, 


y->Kx)  [cp  G  §(x)i^^^^  <  1  V<j>3. 

y'"(x)  the  set  of  virtual  mole  fractions  associated 


There  is  no  reason  to  suppose  that  M(f1h)  consists 
of  but  one  poliit.  Problems  exist  vjhich  possess  no  solu¬ 
tion,  and  others  can  be  constructed  V7hich  possess  many. 
Hov7ever,  whenever  M(F|H)  is  non-empty  and  H(A,  b)  contains  a 
y  >  0,  V7e  can  prove  uniqueness  of  a  sort. 
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Theorem  II. 6:  Suppose  3  z  e  H(A,  b)  3  z  >  0«  Then  for 

any  x,  y  e  M(F|H),  ?*(x)  «  $*(y).  That  is,  there  is  a  unique 

set  §*  of. virtual  mole  fractions  associated  with  problem  (1.6), 

icl 

Proof:  Define  the  carrier  of  a  vector  x  €  E  ,  x  >  0  by: 

C  (x)  “  { j  I X j  ^  0 }  . 

Earlier  in  this  section  vje  defined  a  set  J  of  indices  of 
variables  in  compartments  which  did  not  vanish,  and  in 
the  course  of  Theorem  II. 4  argued  that  if  x  e  M(F[H),  then 
j  e  J  =>  Xj  >  0.  Thus  for  x  e  M(FjH),  J  ~  C(x). 

Lemma  9.7,  pg.  370  of  [5J  states  that  if  M(F1H) 
is  non-empty,  then  there  exists  x°  e  M(FlH)  such  that: 

C(x)  c  C(x°)  y  X  G  M(FlH)  . 

From  this  and  [5]  Lemma  9.5,  pg.  370,  (which  states  that  if 
X,  y  G  M(FiK),  then  x^  =  y.  whenever  both  are  defined)  we 
know  that  for  every  x  g  M(F1H), 

(11.16)  X.  =  x9  V  j  G  C(x)  . 

J 

Clearly,  then, 

(11.17)  «>:  (x°)  c  $-x  (x)  y  X  G  M(FlH)  . 

Suppose  o  G  $>Kx) .  Let  8  =  x°  —  x.  Let  us  now  eval- 

j 

uate  Fq(x).  We  may  do  this  s5.nce  clearly  0^  >  0  if  x^  =  0, 
and  hence  0  is  admissible  at  x.  Further,  x,  x°  e  M(FlH). 


lU 


Thus : 


(n.i8)  f'(x)  -  S  ■  e.-Co.  +  log  ft,) 

“  JeC(x)  J  J  j 

« '• 

+  Z/9j  (Cj  +  log  Sj)  •  0. 
jeC(x®)-C(x) 

On  the  other  hand,  cp  e  $*(x)  tells  us  that: 

(11.19)  D  0g(c.  +  log  ep.)  »  0  . 

jeC(x°)  J  J 

Subtracting  (11.18)  from  (11.19)  we  find'  that; 

.  (11.20) •  0  a  D  0j  (log  cpj  -•  log  0j) 

jeC(x°)-C(x) 

^  Doj  iog  v<j>  <  0. 
jeC(x°)-C(x) 

The  first  inequality  holds  as  an  equality  if  and  only  if 
0j  a  cp^  for  each  j,  by  Lemma  II.  5.  The  second  inequality 
holds  as  an  equality  if  and  only  if  1  each  <j>. 

That  is,  cp  e  §*(x)  'Pj  ^ j  ^  j  e  C(x°) .  Hence 

cp  e  $*(x°)  .  Thus 

(11.21)  $*(x)  c  §*(x°). 

Combine  (11.21)  with  (11.17)  and  the  theorem  is  proved.  Q.E.D. 

We  can  more  fully  describe  $*  as  a. consequence  of  the 
next  several  results.  We  let; 

P*  a  {xT  e  e’'”  I  c  +  log  cp  =3  A^rf  for  some  cp  e  5  ] . 
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Lsmma  11,7! 


■  P*  is  convex. 


?roof  i  Let  x  e  M(F|H(A,  b))  have  a  maximal  carrier  J  =  C(x) . 
Then  we  may  say  that  n  e  P*  if  and  only  if  tr  satisfies! 

!A^jTT  «  Cj  +  log  Xj 

. 

exp(A?TT  -  c. )  <  1  for  each  j  ^  J. 
ke<j>  ^ 

If  we  thei.  let  to.  “  expCATn  -  c.)  for  any  n  satisfying 

J  J  J 

(11.22),  the  vector  tp  e  $*. 

Since  'exp*  is  a  convex  function,  it  is  clear  that 
the  constraints  (11.22)  describe  a  convex  set.  Q.E.D. 

Lemma  II. 8:  Suppose  there  exists  y  e  H(A,  b)  such  that 
y  >  0,  and  let  the  matrix  A  have  full  rank  m.  Then  P* 
is  compact. 

Proof !  P*,  the  set  of  tt  satisfying  (11,22),  is  surely  closed. 

If  it  is  not  compact,  it  must  be  unbounded. 

Since  P->^-  is  convex,  if  it  is  to  be  unbounded,  by  [4], 
Lemma  3,  there  must  exist  a  ray  0  such  that  for  any 
n  e  P*,  n  +  t??  e  P*  for  every  t  >  0.  Clearly  (11,22)  implies 

if  j  e  J 

if  j  J 


that  71  must  satisfy: 


(11.23) 


aT?;  =  0 


A??!  <  0 
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Let  X  e  M(F|H)  as  in  Lemma  II.  7,  and  let  y  e  H(A,  b) 
satisfy  y>0.  If0»y  —  x>  then  clearly  0j  >  0  if  j  ^  J. 
Further,  since  x  and  y  e  H(A,  b),  A0  «  0.  .  Thus,  by  (11.23), 


(11.24) ,  eVn  “  S  0.aT??  «  S  0.aT??  “  0. 

j«l  J  J  J 

Again  by  (11.23),  since  0j  >  0  <  0  if  j  ^  J,  we  must 

have  «  0  if  j  J.  Hence, 

(11.25)  A^??  =  0. 


Since  A  has  full  row  rank  m,  (11.25)  implie^sf'tlwjr  ? 


(11.26)  7?  »  0. 

Thus  P*  is  bounded.  Since  it  is  closed  as  well,  P*  is 
compact,  Q.E.D. 

Corollary  II.  9;  Q*  =  {a  e  |  a  «  A^n  for  some  n  e  ‘P*) 
is  compact  if  there  exists  y  e  H(A,  b)  satisfying  y  >  0. 
That  is,  Q*  is  compact  whatever  the  rank  of  A. 

Proof;  Let  A  have  rank  r  <  m,  and  let  B  be  an  m  x  m  non- 
singular  matrix  such  that  B^A  “  (  q  j '  X  is  an  r  x  n 

matrix  of  rank  r,  and  "0"  denotes  the  (m  —  r)  x  n  matrix 
consisting  entirely  of  zeros.  Let|^-  B^b,  Clearly, 

M(FlH(A,  b))  =  M(FlH(X,  K)). 

Define  P*  to  be  the  set  of  r— vectors  7?  such  that; 


V  j  J. 


XjTJ  “  Cj  +  log  Xj 

S  exp  ^1,)  <  ^ 

ke<j>  ^ 

Clearly,'  e  ?*  if  and  only  if  for  some  arbitrary  (m  ~  r)- 
vector  y,  «  n  e  P*. 

Xi6t^  • 

«  [a  e  E^ja  »  for  some  71  e  F*} . 

By  Lemma  II. 8,  P*  is  compact.  As  the  range  of  a  con¬ 
tinuous  (in  fact  linear)  function  whose  domain  is  compact, 
surely  is  compact.  But  it  is  easy  to  see  that  Q*  =  Q*. 

Q.E.D. 

Theorem  II. 10;  $*  is  compact  if  there  exists  y  e  H(A,  b) 

satisfying  y  >  0, 

Proof ;  CO  e  if  and  only  if  co.  =  exp(a.  —  c.)  for  some 
a  e  Q*.  Since  is  the  range  of  a  continuous  function 
whose  domain  Qx-  is  compact,  $*  must  itself  bo  compact. 

Q.E.D. 


It  is  not  true  that  is  a  convex  set  in  general. 

We  defined  §*  in  such  a  way  that  for  any  compartment 

<j>  which  need  not  be  empty  at  equilibriam  (i.e.,  such 

that  for  some  x  e  M(F|H),  >0),  1  for  every 

1  2 

?  e  $*,  and  further,  any  §  ,  ?  e  agree  in  compartment 
<j>.  The  next  result  extends  this  idea  slightly.  It  will 
prove  important  later. 
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■Lemma  11,11:  Suppose  for  every  §  e  $*,  *»  1  for  a  par¬ 
ticular  <j>*  IThen  for  eacfe  k  e  <j>,  and  every  e  $*, 

■pl  = 

^k  V 

t  *9  19 

Proof ;  Suppose  §  ,  §  e  $*4  and  for  some  k  e  <j>. 

Let  ■TT^,  e  P*,  be  such  that: 


exp(A^Tr^  —  c) 


i  «  2, 


Since  P*  is  convex  (Lemma 
of  P*.  Let: 


II.  7),  ??  «  H — is  an  element 


Y  «  exp  -  c)  . 


Clearly,  then. 


kj> 


1.2 1  * 


ke<j> 


k'^kj 


But  for  any  a,  b  >  0,  a  5^  b,  it  is  well  known  that: 

(ab)  ^  < 


^  a  +  b 


Thus 


1  I  k 

,  since  for  some  k  c  <j>,  5^  f  §2^  have: 


X 

s 


'f^  4.  p2' 

^k  +  ^k 


=  s  (§^4)  <  ^  s  1 


'<J>  ke<j> 


ke<j> 


QED. 


IQ 


HI;  Overcoming  Deaeneracy  With  Slacks 

Although  we  now  ’icnow  how  to  recognize  whether  a  com¬ 
position  X  is  a  solution  to  problem  (1.6)  we  still  face 
several  difficulties  in  trying  to  solve  the  problem.  One 
of  them  is  constructing  the  vector  which  satisfies  the 
assumption  made  throughout  the  previous  section^  that  there 
exists  a  vector  x  e  H(A,  b)  such  that  x  >  0.  This  is  discussed 
in  [7]. 

Other  difficulties  concern  degeneracy.  It  may  be, 
for  example,  that  the  solution  set  M(FlK'  of  a  particu].ar 
problem  contains  more  than  one  vector  x.  Or  it  may  be  that 
for  some  x  e  M(F[H),  at  least  one  phase  i;.  c-inty.  To  cope 
with  this  last  difficulty,  v;e  have  construct'  d  the  <.  Laborate 
theory  of  Section  II.  This  chapter  develops  from  that  theory  a 
method  that  avoids  both  of  these  difficulties. 

1.  Slacks;  General  Fo 

Ue  may  insure  in  no;  •.  cases  that  the  problem  we  ac¬ 
tually  bolvo  ’..’ill  have  a  unique,  strictly  positive  solution 
by  addin y  on^  extra  variable  to  each  compartment.  Those 
extra  variables  are  called  slacks .  With  them,  the  problem 
becaries : 


X. 


+  log  ~ 


X 


+  S, 


<j>  "  xj> 


+  E 

<j> 


log - - \ 

\3>  ^<3>l 


(III. 1.1) 


s.t.  Ax  »  b 


>  0 


X  ^  0 


The  mole  fractions  are  of  course,  changed  by  the  addition 
of  slacks.  Thus  in  (III. 1.1) 


(III.  1.2)  =  r 


X. 


\3>  ^<j> 


<j> 


X  -  +  S  . .. 

<J>  <J> 


where  is  as  it  was  before. 

Notice  that  this  method  of  handling  degeneracy  does 
not  require  that  the  problem  actually  be  expanded.  Instead, 
we  may  treat  the  slacks  implicitly,  choosing  a  small  pos¬ 
itive  for  each  compartment  <k>  and  solving: 

"<j>  k3>) 

(111.1,3)  s.t.  Ax  =  b 

X  ^  0 

A 

To  compute  and  in  (III. 1.3)  we  have  substituted 

^<j>  ^<j>  latter  occured  in  (111.1.2). 

We  define: 
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(III.  1.4)  Fg(x) 


=  S  x.(c.  +  log  X.)  +  E 

j  J  J  J  <j> 


"<j>  ®<j> 


xdiere  fi  is  the  vector  whose  components  are  the  fi 


'<j>' 


Theorem  III, 1.1;  Assume  there  exists  y  e  H(A,  b)  satisfy¬ 
ing  y  >  0.  For  each  £  >  0,  let  M(Fg]H(A,  b))  be  nonempty. 
Then  it  contains  exactly  one  point  x(e),  and  that  point  is 
a  strictly  positive  vector. 


rroof ;  Suppose  M(F^lH)  is  non-empty.  Looking  to  (III. 1.1) 
we  see  that  by  adding  slacks  to  the  problem,  we  have  in¬ 
sured  that  the  sum  of  variables,  including  the  slacks,  for 
each  phase,  must  be  strictly  positive  (in  fact  it  can  be 
no  smaller  than  >  0) .  Thus  the  set  J  (of  indices  in 

non-vanishing  compartments)  must  always  include  every  index 
# 

J* 

But  v:e  know  frar-  Theorem  II. 4  that  the  solution  must 
be  strictly  positive  for  each  index  in  J,  Hence  every 
X  e  M(Fg|H)  is  strictly  positive. 

Lemma  9.5,  p.  370,  of  [5]  states  that  if 
•X,  y  6  M(F^'/.)^  then  since  x,  y  >  0: 


"  ^3"  j  ®  2,  ...,  n 

V7here  x^  and  y^  arc  defined  according  to  equation  (III. 1.2). 
In  particular: 


\3>  ^<j> 


:<3> 


^<j>  ^<j> 
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Thus 

X  «  y.  QED. 


Simple  algebra  shows  from  this  that 


Theorem  111,1.2;  Assume  there  exists  y  e  H(A,  b)  such  that 
y  >  0,  For  each  6  >  0,  M(Fg|H)  is  non-empty  if  and  only 
if  M(F|H)  (the  solution  set  .of  the  problem  unperturbed  by 
slacks)  is  bounded  and  non-empty. 

Proof:  (=»)  Suppose  M(F|H)  is  empty  or  unbounded.  Then 
[4].,  Theorems  2^  4  and  5,  state  that  there  must  be  a  non¬ 
zero  vector  0  e  satisfying: 

F(e)  ^  0 

(III. 1. 5)  A0  =  0 

0  >  0,  and  0  'j’  0. 

Choose  any  y  e  H(A,  b)  such  that  y  >  0.  (We  know  from 
the  previous  Theorem  III. 1.1  that  only  such  y  are  candidates 
for  a  solution  to  the  perturbed  problems.)  Then: 

F(y  +  t0)  <  F(y)  +  tF(0) 

by  the  convexity  and  homogeneity  of  F.  Thus 

(in.1.6)  F'(y)  =  Urn,  ^  P(0)  <  0  . 

”  t-'O 


From  equation  (11,3)  we  may  compute  the  directional 
derivative  of  F^  at  y  in  the  direction  0  to  be: 

(III,  1.7)  ^  —  . 

t"'0  <i'>  ''J'  V  .  4-  p 

^<j>  ^<j> 


2^ 


Since  e  >  0,  each  of  the  sums  nonnegative.  Because 

0  0,  for  at  least  one  <j>  we  must  have  >  0«  For 

every  the  fact  that  >  0  insures  that: 

(ni.1.8)  log  — —  <  log  1  =  0. 

y<j>  ®<j> 

Combining  (VI. 1.6-8)  we  find  that: 

(III.  1.9)  llm^  ^  F^(y  +  te)  <0. 

Thus  by  Theorem  II. 2,  y  4  M(Fg|H).  Since  y  >  0  was  arbitrary 
and  since  y  e  M(Fg|H)  implies  y  >  0,  we  Icnow  M(Fg|H)  is  empty 

(<=)  If  M(F^|H)  is  empty,  then  by  [4],  Theorems  2,  4  and  5, 
there  must  exist  a  vector  0  satisfying  (III. 1.5).  This 
in  turn  implies,  by  [4]  Theorems  4  and  5,  that  M(F|H)  is 
either  unbounded  or  empty.  Q.E.D. 

We  V7ill  assume  throughout  the  remainder  of  this  chapter 
that  M(F|H)  is  bounded  and  nonempty.  In  addition  we  will 
suppose  that  there  exists  y  e  H(A,  b)  satisfying  y  >  0. 

This  will  save  us  from  including  these  assumptions  in  the 
hypotheses  of  each  theorem. 

Theorem  III.  1.3;  For  each  vector  fi  >  0,  let  M(F^|II(A,  b)) 

=  f2i(c)}*  Then  x(p.)  and  x(r.)  are  continuous  functions  of 

e  >  0. 

Note:  This  statement  is  equivalent  to  saying  that  x(f,) 
solves  (III.  1,1).  By  the  previous  tv70  theorems,  x(f)  is  a 
positive,  single-valued  vector  function  of  the  vector  B. 


Proof;  Since  x(e)  >  0,  x(e)  is  well  defined  and  continuous 


if  x(6)  is  continuous.  To  see  that  x(G)  is  continuous, 
see  [8],  Cor.  1.2.4,  p.  528.  QED, 

It  is  our  purpose  in  this  section  to  show  that  the 
continuity  of  x(6)  and  x(G)  extend  in  some  sense  to  6  =  0, 
which  characterizes  the  original  problem,  (1.6).  That 
is,  we  wish  to  show  that  if  each  component  of  G  is  chosen 
sufficiently  small,  then  the  solution  x(G)  will  be  as  close 
as  desired  to  some  solution  of  the  original  problem,  and 
x(G)  will  approximate  a  vector  of  virtual  mole  fractions. 

Lemma  III,  1,4:  Let  {g^^3  be  a  sequence  satisfying  G^  >  0  vk 

and  G^'"  -  0.  Then  there  exists  an  infinite  subsequence 

S  of  (i,  2,  ...)  such  that  x°  =  lim  x(C^^)  exists,  and  for 

keS 

each  such  subsequence, 

(Ill.l.iO)  x°  G  M(FiH(A,  b)). 

Proof:  [83,  Corollary  II. 3.1,  p.  545,  and  Theorem  1.2.2, 

p.  526;  and  Theorems  III. 1.1-2  above.  Q.E.D, 

Thaore-m  III.  1.5:  Let  {G^^}  be  as  in  Lemma  III. 1.4. 
exists  a  subsequence  S  of  (1,  2,  ...)  such  that  % 
exists,  and  for  each  such  §  we  have  that: 


Then  there 

=  lim  >c(6'^) 
keS 


(III. 1.11) 


Proof:  Since  for  every  j  and  k,  0  <  x.(e^)  <  1,  we  can 

^  V 

find  a  subsequence  S  of  (1,  2,  ...)  so  that  §  «  lim  x(0 

keS 

exists. 

By  Lemma  III.  1.4,  we  may  find  a  subsequence  S’  of  S 

such  that  x°  =  lim  x(e^^)  exists,  and  x°  e  M(F|H). 
lC€  S  * 

Since  x(e^)  e  M(Fg|H),  and  since  x(e^)  >  0  for  each 
k  (making  every  reaction  vector  admissible  at  x(e^)),  we 
know  that  for  each  0  satisfying  A0  =  0,  we  must  have 
0»(c  +  log  x(e^))  =  0.  Taking  the  limit  of  this  expression 
as  k  €  S’  becomes  infinite,  we  find: 

(m.1.12)  0.(c  +  log  §)  =  0  V0  3  A0  =  0. 

Clearly,  for  each  phase  <j>, 

(ni.1.13)  T,.  =  Hm  (  E  xfe*'))  <  1. 

keS'  te<j> 

It  is  easy  to  show  that  for  each  i  such  that  X?  is 

J 

defined. 


(III. 1.14) 


e  ^  O 

§ .  «  X. 

J  J 


and  that  for  every  j. 


(III. 1.15) 


5j>0. 


By  (111.1.12—15),  and  Theorem  II. 6, 


§  €  ^  . 


Q.E.D. 


2,  Slacks;  Special  Form 

The  mp.st  natural  sequences  one  may  consider  a^^e 

those  which  take  the  form, 

(III. 2.1)  “  ty.c 


where  a  is  a  constant  veccor  with  one  strictly  positive 


element  corresponding  to  each  compartment  and  is  a 

sequence  of  positive  real  numbers  whose  limit  is  zero. 

In  this  section  V7e  will  explore  the  properties  of  the 
sequence  of  solutions  to  problem  (III. 1.3),  where 

takes  the  form  of  (III. 2.1). 

We  shall  change  our  notation  somewhat.  Instead  of 
writing  Fg(x),  we  will  \^ite  F(xj  t),  where: 

r  /  X.  \ 


(III.2.2)  F(x;  t)  - 


E  X, 


c.  +  log 


JL 


^<j> 


« 


to 


+  Z 

<j> 


ta 


<j> 


log  - 


j>  J 


The  sums  x. ..  do  not  include  the  slacks,  ta^..  .  Also, 

<3> 

in  witing  F(x;  t)  x-ze  understand  that  the  vector  o  will 
be  constant,  and  so  do  not  mention  it.  In  the  same  x^/ay, 
x^e  Xi7ill  replace  x(ta)  by  x(t),  leaving  the  constant  vector 
a  out. 

We  shall  show  that  x(t)  approaches  a  limit  as  t  -♦  0, 
and  Xi7e  shall  identify  the  limit.  We  shall  attempt  to  do 
the  same  for  x(t),  but,  as  the  reader  x^ill  see,  x-je  fall 
somewhat  short  of  complete  success. 


0*7 


Take  x(t)  first.  Define  the  sets: 


J  «  e  5*3 

I  =  {j  i  J}. 

In  addition  it  will  be  convenient  to  define  the  related 
sets: 


J*  “  {<j>lj  e  vT)  =  ^  s  5*3 

I*  =  (<j>jj  G  13  =  {<j>  i  J*3 


We  know  that  for  j  e  J,  Xj(t)  must  converge  as  t  0, 

since  it  is  bounded^  since  by  Theoran  111,1,5  every  limit 

point  of  x(t)  must  be  an  element  of  and  since  by  Lemma 

II. 11,  there  is  for  j  e  J  only  one  possible  value  for  §. 

J 

if  §  €  §*. 

For  each  j,  define  a  function  of  tt  as  follows: 


(III. 2.3)  5:j(tr)  =  expCA^TT  -  c^  . 

We  know  from  Theorem  ,11.4,  and  definitions  (II. 5b)  and  (11.15) 
that  if  TT(t)  is  the  vector  of  Lagrange  multipl5.ers  cor¬ 
responding  to  the  solution  x(t),  then  x.(t)  =  §.(n(t)). 

*3  J 

Substituting  for  Xj(t)  from  equation  (III. 1.2),  we  may  solve 
for  the  quantities  Xj(t),  1  <  j  f  ^  in  terms  of  the  quantities 
5:i(TT(t)),  l<j<n  and  the  slacks  ta  •  The  result  is 
just: 


(III. 2.4) 


t:a^^^§:j(rr(t)) 
~T<j^(n(t))  ’ 


?8 


Suppose  that  e  M(F|H).  Then  we  Imow  that: 


(III.2.5)  2  A.x.(t)=~S  A.(x.(t)  -  X?). 

jel  -3  J  jeJ  J  ^  ^ 


Substituting  (III. 2.4)  into  (III.2.,5)  for  j  e  I  we  find  that: 


tcr^^v.-^4§^(Tr(t)) 

(III.2.6)  E  — ■  =  _ 


jel  l  -T<j^(,T(t)) 


1  A.(x.  (t)  —  X.) » 

jsJ  ^  ^  J 


A  simple  computation  from  (III. 2.3)  shows  that  if: 
g^(TT)  =  ~  S  _  log  (1  ~  T^^^(tt)) 


<j>el* 


<j> 


then: 
(III. 2. 7) 


ta^.^A.§.(!T) 

S  - -  «  vg  (rr)  , 

jel  (1  -  (rr)) 


By  (III. 2. 7),,  equatioii  (III. 2. 6)  is  the  optimality  condi¬ 
tion  for  the  problem: 


(HI.  2, 8) 


Min  g^(-,T) 

T 

S.  t.  AjTT  =  Cj  +  log  Xj(t)  , 


That  iSj  if  we  knew  for  any  fixed  t  the  values  of  Xy(t), 

*J 

then  T7(t)  would  be  a  solution  of  (III. 2. 8). 


Lemma  III. 2.1:  If  the  matrix  A  has  full  rank  m,  then  (III. 2. 8) 
possesses  a  unique  solution  rr(t) . 

Proof:  We  know  that  (111,2.8)  is  feasible.  That  is  rT(t) 

T 

satisfies  AjTi(t)  =  Cj.  +  log  Xj(t)  . 
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For  t  >  0,  we  niey  compute  the  matrix  of  second  partial 


derivatives  of  g. 
(III. 2. 9) 


This  becomes, 

j€l  (1  -  <j>el*  (1  - 


\  . 

where  L.  “  S  A,  Clearly  v^g  (n)  is  a  positive 
.ke<j>  ^ 

semidefinite  matrix  for  all  tt  such  that  j  ^  I, 

so  that  g(TT)  is  convex  on  P*. 

T^tn 

Further,  for  any  v  e  E  , 

(III. 2. 10)  v'^(v2g^(Tr))v  =  0  «  AjV  =  0. 


To  see  that  (III. 2. 10)  is  true,  let  Y  be  a  diagonal 
matrix  with  as  many  columns  as  there  are  indices  in  I.  The 


jj^-1  component  of  Y  shall  be  y^  =  ‘  ^ 

^  ■"  ^<j> 


>  0.  Then  the 


first  term  of  (III. 2. 9)  becomes  A^YA^;  by  elementary  linear 

mm  .r.  aT 

algebra  one  may  show  that  v'^A^^A^v  -  0  if  and  only  if  A^v 
A  similar  comnient  applies  to  the  second  term  of  (III. 2. 9)  . 


1  9 

If  there  are  two  solutions  n  and  to  (III. 2. 8)  they  must 
differ  by  such  a  v.  Further,  since  each  solution  satisfies 
the  constraints,  -  rr^  =  v  must  also  satisfy  AjV  =  0. 

Hence  Av  =  0.  Since  A  has  full  rank  m,  v  =  0,  so  the  solu¬ 
tion  n(t)  of  (III. 2. 8)  is  unique.  Q.E.D. 


Next  V7c  wish  to  establish  that  x(t)  must  approach  a 
limit  as  t  -♦  0.  We  v?ill  again  assume  that  A  has  full  rank 

m. 
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Theorem  111,2.2;  §  lim  x(t)  exists,  and  §  «=  ?(n*),  where 

t-^0 

IT*  is  the  unique  solution  to: 


(III. 2. 11)  Min 


-  S  CT...  log  (1  -  T<4>(TT))j 

L  J 


m 

S.t,  AjTT  =  Cj  +  log  ?j  . 

(5^  =  iim  x,(t),  which  we  know  exists.) 
t-^0 

In  particular  ^  each  <j>  e  I*. 


Proof :  For  each  t  >  0,  tt  (t)  must  satisfy; 


(III. 2. 12)  Min 


-  S  a  .  log  (1  ~ 
<j>Gl* 


m 

s.t.  AjTT  =  Cj  +  log  x;j(t)  . 


This  is  just  (III, 2. 8)  V7here  we  have  divided  the  objective 
function  g^(TT)  by  t. 

We  l<now  that  2^  (t)  is  continuous  for  j  e  J,  and  is 
bounded  away  from  zero.  Thus  the  right-hand  side 
(Cj  4-  log  >l:j(t))  of  the  constraint  set  of  (III. 2. 12)  is 
continuous  at  t  =0. 

Furthermore,  at  t  =  0,  problem  (III. 2. 12)  has  a  uniqvie 

solution  by  Lemma  III. 2.1.  Thus  by  [8],  Corollary  II. 3.1 

and  Theorem  1.3.2,  vjo  knox7  that  the  solution  Tr(t)  to 

(111,2.12)  is  continuous  at  t  =  0;  hence  x(t)  is  continuous 

at  zero,  and  §  =  lim  x(t)  must  satisfy  (III. 2. 11).  Q.E.D, 

t-0 

Next  wc  wish  to  v7ork  on  :c(t)  itself.  As  wc  stated 
earlier,  v/c  have  l)con  unable  to  prove  that  x(t)  approaches 


a  limit  as  t  -*  0.  For  the  reasons  that  folio  >7,  we  conjecture 
that  x(t)  does  have  a  limit. 

Let  §  e  ?*  be  the  limit  of  x(t)  as  described  in  Theorem 
III. 2. 2.  By  [5]j  Lemma  9.7,  p.  370,  we  Icnov;  there  exists  • 
y°  €  M(F|H)  such  that  for  every  compartment  <j>,  “  0 

implies  that  “  0  each  x  e  M(F|H),  We  define  the 

index  sets: 

i  i*=C<j>lj6i) 

J  =  tj|y<j>  >  j  j*  =  {<j>|j  s  j) 

K  =  (j  4  J|5<j^  =  1)  ;  K*  =  t<j>lj  e  K)  . 

Lemma  III. 2. 3:  Let  M(F|H)  be  bounded  and  nonempty,  and 
suppose  there  exists  y  e  H(A,  b)  such  that  y  >  0.  Then 
there  exists  a  unique  solution  x*  to  the  problem: 

'(111.2.13)  Min  gOO  =  Min/-  S  o...  log 

v  <j>6 


s.t,  X  p  M(FlH)  . 


Further,  for  each  j  e  J,  x*j  >  0, 

Proof:  From  [5],  Theorem  9.5,  p.  370,  we  can  write  M(F|H(A,  b)) 


t  s  the  set  of  solutions  to  a  system  of  linear  equations  together 
with  nonnegativity  constraints.  Let  y°  e  M(F|H)  satisfy 
yj  >  0  for  each  j  e  J,  and  define: 


(III. 2. 14) 


E 

lCG<j> 


j  e  J. 


Then  x  g  M(F|H)  if  and  only  if: 
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(III. 2. 15) 


(a) 

(b) 

(c) 


Xj  =  yj  for  every  j  e  such  that  >  0^ 

X.  =  0  for  j  ^  J 
J  ^ 

<j>^  J<j>^<j>  “ ^<j>  ^  “• 


If  we  find  any  solution  x  to  (III. 2. 15c)  we  may  reconstruct 
a  unique  element  x  of  M(F|H(A,  b))  by  using  (III.2,15a--b)  . 
Further,  that  element  x  V7ill  satisfy  x ^  >  0  for  each  j  e  J 
if  and  only  if  >  0  for  every  j  e  J. 

Thus  (III. 2. 13)  is  equivalent  to  solving: 

(III. 2. 16)  Min  g(x)  =  Min  /-E  a.,  log  BL-v] 

\<j>eJ-x-<J>  <J>/ 


s.t. 


S 

<j>GJ 


.?<3>==<j> 


=  b 


It  is  easy  to  see  that  g(x),  as  a  function  of  the  x’s, 
is  strictly  convex.  Since  the  constraint  set  is  bounded 
(i.e.,  M(FjK)  is  bounded),  (III. 2. 16)  has  Ji  vniijre  solu¬ 
tion.  Since  there  exists  a  positive  feasible  solution, 
and  since  g(x)  is  infinite  if  some  =  0,  j  £  J,  It 

is  clear  that  the  solution  x'"  to  (III. 2.16)  is  positive  on 
J.  Q.F.D. 


We  conjecture  that  x(t)  has  a  limit  as  t  -  0,  and 

that  lim  x(t)  =  x*. 
t-^0 


•  By  definition,  x(t)  is  the  unique  ■  oint  in  M(F(  •  ,  t)lH). 
Thus,  x(t)  is  also  a  solution  (though  not  necessarily  the 
unique  solution)  to; 

(III. 2. 17) 


Min  G(y,  t)  =  Min  B  c,  +  log  — ^  -f  log  ~ 

lUK  J\  J  .  Y...  X 


y<j> 


+  B  y.  c .  +  log  ~ 

T  J  \  J  V 


y<j>  +  ‘’<j> 


s,  t.  ye  H(A,  b)  . 

It  is  easy  to  see  that  G(y,  t)  is  convex  for  fixed  t,  though 
it  is  not  necessarily  strictly  coin/ex;  and  surely  x(t)  sat¬ 
isfies  the  condition  that: 

»  o.  +  log  x,(t)  =  ATn(t). 

J  J  J 

We  separate  G(y,  t)  into  two  parts.  Let: 


(ITI.  2.13a)  g(y,  t,  p)  =  r  +  log  ^ 

'  ^<j>‘ 


X 


+  PT:  y  .  log  - - 

<j>elx-  +  tcr^.^ 


(III.2.18b)  h(y,  t,  p)  =  i[G(y,  t)  -  g(y,  t,  p)]. 


If  we  then  define  the  set: 


(III.2.19)  Sp(t)  =  {y  e  Il|g(y,  t,  p)  ^  g(x(t),  t,  P) } 


then  it  is  obvious  that  for  any  t  >  0, 


(11112.20)  y  €  Sp(t)  =»  h(y,  t,  p)  ^  h(x(t) .  t,  p) . 
Equivalently  we  may  vnrite: 

(111.2.21)  x(t)  e  M(h(-,  t,  ^)jS^(t):), 

Leipma  III. 2.4:  Let  M(F(*,  0)  }H)  be  bounded  and  nonempty.- 
Then  for  any  0  <  p  <  1,  M(F(*,  0)  |H)  c  S^Ct) ^  for  all  t  >  0 
sufficiently  small. 


Proof:  Note  that  g(y_,  tj  p)  is  convex  in  y  for  fixed  t 
aiid  p.  Let  x°  e  M(F(*,  0)  l.H). 

Notice  that  ^(y,  tj  p)  is  a  G5.bbs  function  with  each 


Cj,  for  j  £  modified  to: 

t 

^  _  _ _ ... 


?j  “  Cj  +  p  log 


-f  p  log  i  e  j: 


Let  §  •-  lim  x(t);.  as  in  Theorem  111.2.2.  We  know  that 
t-‘0 

5  e  f*.  If  V7e  replace  c.  by  d.  for  every  j  g  I,  where; 

•J  mJ 


d .  ==  c  .  log  k  ^ .  i  G  I 

3  1  <J>  ** 

and  log  >  log  then  clearly  §*  defined  by 

is  an  element  of  the  new  so  that  x°  is  a  solution  to 
the-  perturbed  problem. 

But  for  given  p  <  V7e  can  alv7ays  pick 
Cg  >  0  sufficiently  small  so  that  for  0  <  t  <  t^, 

W  ~  p' 
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P  log 


_ >  log 


Thus  X®  €  Sp(t)  for  0  <  t  <  tp.  Since  g(x,  t  P)  is  cbiv- 
stant  oh  M(F(«,  0)1  H), 

M(F(*,  0)Th)  c  Sp(t)  0  <  t  <  tp. 

Q.E.D. 

The  next  item  we  must  study  is  the  function  h(y,  t,  p) , 
This  we  deal  with  in  two  leinmas. 

Lemma  III ..2. 5:  Let  x®  €  M(F(»,  0)jH)  be  fixed.  Then 


-E  a...  log  -2  a<.  if  C(x°) 

SJ-^  ^4**e.T* 


<j:>eJ* 


lim  h(x°,  t,  P)  « 
t-‘0 


if  C(x°)  c  J. 


Proof;  Since  x°  e  M(F(»,  0)lH),  x?  «  0  if  j  e  I  U  K. 
Thus  from  (III. 2. 17)  and  (III. 2. 18b), 


h(x°,  t,  «  ,  4  E  log  - 


<j>  “’<j> 


Suppose  for  some  j  e  J,  Then: 

x°  log  - - -  «  0  by  definition,  for  all  t  >  0 

j. 

X  +  to^j^ 
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. 


and: 


-^<i>  "<j>  “  * "  ‘ 


On  the  other  hand,  if  >  0#  we  see  that: 

-o  T  “^-o  /_  \k 


1 108 


Also, 


^<3>  1-  ^ 

*°<J>  ■*■  *”<J>J  V  <j)/ 


<j>‘ 


-o<j>  log  (3c»^j>  +  to^j>)  -  -a<j>  log 
Thus  the  lensna  is  proved.  QED. 

Lemma  III, 2. 6:  Let  {t^)  be  a  sequence  satisfying  >  0  ^ 

for  every  n,  “*  0  and  x(tj^)  -•  y°  €  M(F|H)  as  n  -*  «. 

Let  h^  «  h(x(tj^) ,  t^,  p) ..  Then  if  the  carrier  C(y®)  ^  J, 

lim  h  ».  Otherwise  (i.e.,  if  C(y°)  «  J) 
n-*»  “ 

I 

lim  h„  «  (1  -  P)  2  J<j>  «  E  a ... 


1  -  T. 


"<j> 

'^j>6K*UJ* 


-  2  a  ,  log  y  . 
where  5  lim  x(0 . 

n-*eo  “ 

Proof:  From  equation  (111.2.4-),  if  j  e  I, 


(III.2.22)  11m  k.,.  (t)  log  -i 

nwco  ^  NJ^  V' 


1  -T, 
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If  j  c  cJ  U  kv  We  see  that: 


log  — 


to 


+  0 


to 


X^4v  +  to^av  Xviv  +  ta 


Xj> 


<j>  Xj> 


<j> 


^*<j>  +  '=’<jy 


X.iv  (t  ) 

Multiplying  by  and  taking  the  limit  as  n  -♦  », 

^n 

clearly  the  error  term  vanishes.  We  are  left  with: 

(HI.  2. 23)  lim  log - - — 

n^»  n 

-  ( j  €  J  U  K) . 


The  remaining  term  is  easy. 


(III.2.24)  -lim  E  a...  log  (x...(tj  +  ta  a) 
n-*»  <j>€ 


If  C(-y®y  ...J 

-  if  C(y°)  c  j. 

Putting  Equations  (111.2.22-24)  and  (III. 2. 18b)  together, 
yield  the  result  we  wish,  QED. 

These  last  three  lemmas  enable  us  to  prove  the  following 
theorem. 


Theorem  HI. 2. 7:  Let  M(F(*,  0)|H)  be  bounded;  Let  {t^^}  be 
such  that  t^  >  0  for  each  n,  t^^  0  as  n  -♦  «,  and  x(t^) 
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-■y°  e  0)lH)  as  n  -  Then: 


01-S 


<j>6j*^J^  <j>CK«r 


Proof;  For  sufficiently  large  n  (i.e,,  sufficiently  small 
t  >  0),  Lemma  III, 2.4  tells  us  that 

X*  G  Sp(t) 


since  x*  e  M(F(«,  0)jH).  Thus  for  n  sufficiently  large^ 
we  have  by  (III. 2.20)  that: 


(III.2.25)  h(x*,  t^,  p)  >  h(x(y,  t^,  p). 


Thus,  taking  limits  as  n  -♦  «  of  both  sides  of  (III. 2. 25) 
we  find  that  Lemmas  VI.2.5-6  imply; 


(III.2.26)  ~  L  a.,  log  x*  >  (1  -  P)  S  „< .  .I. - ^ 

<j>€  ^  J 


-  2  a  .  -  2  a  .  log  . 

<j>GK*  <j>eJ* 


1  -  T 


<J> 


Clearly,  since  the  left-hand  side  of  (III. 2. 26)  is  finite, 
the  right-hand  side  cannot  be  «;  hence,  C(y°)  »  J,  by  Lemma 
III. 2. 6. 

Further,  by  Lemma  111.2.4,  Equation  (III.2.26),  is  true 
for  every  p  <  1.  Thus  letting  fi  «  (1  -  p)  >  0  be  as  small 
as  we  wish,  we  see  that,  as  the  theorem  says. 


r 


Kin.2.27)  Q^  -t  •  log  S  a,,, 

■*<j> 


The  left-hand  inequality  is  a.  consequence  of  Lerrana  XXX.  2. 3^ 
since  y®  .€.M(F'{*,  0)  jH).  QED. 

The  Important  corollaries  of  Theorem  XIX. 2, 7  are  as 
follows . 


Corollary  XXI. 2. 8;  If  y°  e  M(F(»,  0)  jH)  is  a  limit  point  of 
x(t)  as  t  -♦  0,  then  C(y^)  *  J. 


Proof:  Immediate, 


Corollary  XII. 2. 9;  If  K  is  empty  (i.e,,  if  for  every.  ,j  ^  J, 
<  1)/  then  x(t)  has  a  limit  x*  as  t  -*  0. 

Proof;  If  K  is  empty  then  the  upper  limit  on 

^°<i> 

“2  O  ,  log 

is  zero.  Thus: 


0 

<j>e J* 


<j> 


log 


y° 


<i>  “  5*<J> 


By  Lemma  XII. 2. 3,  since  x-^^  uniquely  solves  (III. 2. 13)' 

we  must  have  y°  «  x*,  for  each  limit  point  y®  of  x(t) . 
-QED. 


ho 


‘ — ' 


We  conjecture  that  regardless  of  vdiether  K  is  empty^ 

x(t)  has  a  limit,  and  lim  x(t)  The  reasons  support- 

t-^ 

ing  this,  conjecture  are  three.  First,  every  example  we 
have  tried  leads  us  to  the  limit  x*.  Second,  we  have  the 
following  result. 

Lemma  111,2,10;  Let  every  compartment  <j>  satisfying  j  6  J 

be  identical  to  every  other.  That  is,  if  jj^  and  j^  are 

any  two  indices  in  J,  with  <jj^>  ^  <52>*  every 

^1  ^  there  is  a  k2  6  <j2>  such  that  columns  and 

k2  of  the  matrix  A  are  identical,  and  c.  *  c.  .  In- this 

^  K2 

case  x(t)  -♦  x*. 

Proof ;  Clearly, 

(III.  2. 28)  c.  +  log  X.  (t)  -  aJ  n(t)  “  +  log  (t) 

-  i*  .1  2 

at  the  equilibrium  solution,  since  A.  «  A,  .  Thus:  * 

1  ^2 

(111.2.29)  (t)  B  (t)  , 

^  2 

Summing  over  all  kj^  e  <jj^>  and  e  <j2>,  and  perform¬ 
ing  some  simple  manipulations,  we. find  that  there  exists 
^  pyn^ber  ^v,t)  for  each  t  >  0  such  that: 


(III.2.30) 


+  w<j> 


Kt) 


V  <j>  e  J*. 


If  we  let: 


Ui 


>V 


(III. 2. 31) 


k€<j> 


b(t)  **  b  -  E  A^.Xj^(t)  . 


kelDK 


we  noti<:^;  that  by  (111.2.29),  each  &<jj^(t),,  for  <o>  €  J*, 
is  identical  to  every  otter,  and  that  .since: 


XIII.2.32)  E  P<..(t)(x<..  (t)  +  ta  .  V  »  b(t) 


we  see  that  b(t)  must  be  a -multiple  of  ^^.iv,(t).  Let 

b(t)  =  k(t)p^j^(t).  then  Equation  (111.2.32)  is  equivalent 

to: 


•(III.2.33)  E  (x<..  (t)  +  ta  .  )  =  k(t) 

<j>ej-x-  <3> 


Equations  (III--2i'30.)  say  that  x(t)  solves  the  problem; 


Hin 


(III.  2.34)'  s.t.  E  (x^j;.  +  =  k(t) 


X  e  H(A,  b)  .> 

Since  (III. 2. 32)  and  (111,2.33)  are  equivalent,  problem 
(111.2,34)  is  the  same  as  minimizing  the  same  function  sub¬ 
ject  to  X  G  H(A,  b)  and  x  satisfying  (III. 2. 32)  in  place 
of  (III. 2. 33). 


i»2 


But  x-j(t)  Is  continuous  at  t  «  0,  so  t'hat 
tcohtinuous  at  t  «  0,  and  Xj,jj^(t)  tend  to  zero  a't  it  -  'O,, 
.'SO  that  b(t)  is ■  continuous  and  tends  to  b  at  t  «  0.  Thus 
:in  the  lljnit  as.  t  -  0,  problem  (111.2.34)  is  equivalent  to 
the  problem  that  defines  x-^.,  i.e.,  problem  (HI.  2, 16)  .  Thus 
by  .[8],  Corollary  II.'3..1  and  Theorem  1.3.2, 


lim  x(t)^  «  X*. 
t-0 


QEa). 


The  third  bit  of  evidence  supporting  the  conjecture 
is  that  in  fwo  senses,  a  problem  with  a  nonr-empty  set  of 

indices  K  is  the  limit  of  a  sequence  of  problems  with  K 
empty. 


-Lemma  III. 2. 11;  Let 


jeK 


where  §  e  §*.  .(ByLemmaIl.il,  every?  e  is  the  same 
for  indices  j  e  K.)  Let: 


Hg  <=  {x|Ax  =  b  +  X  >  0). 


Then 


lim^  M(F|Hg)  «  M(F|H). 

NoteJ.:  The  limit  of  a  set  is  as  defined  in  [8], 

Note_2:  The  effect  of  this  perturbation  is  to  move  each 
index  originally  in  K  to  J. 


P;roQf ;  Notice  that  this  perturbation  aliov7s  us  to  con-^ 
struct  a  soiution  x  to  the  problem  from  any  solution 
x°  e  M(F,ill(A,  b))  by 


X.  <= 

J 


(  X?  if  j  6  J 
(  0  if  j  €  I 

(^6§ .  if  j  €  K 


By-  05] j  Theorem  9.5,  p.  370  we  know  that  any  element  of 
M(F|Hg)  must  be  quasidependent  on  this  solution  x.  That 
is,  we  have  for'  the  perturbed  problem  the  same  $*  as  for 
f} i  original. 

But  we  argued  in  Lemma  111,2,3  that  we  could  write  a 
set  of  linear  constraints  describing  M(FIH(A,  b)) ;  and 
we  may  do  the  same  for  M(FlHg.)\  Thus  x  e  M(FlHg)  if  and 
only  if  the  sums  j  e  J  U  K)  satisfy, 

(111,2.35)  S  ^ 


^■<j^  >0  ^ 


where 


ku<j> 


But  by  [8],  Corollary  11.3,1,  the  limit  of  the  set 
of  .X  satisfying  (111,2.35)  is  exactly  the  set  satisfying: 


(III. 2.36) 


<j:>eJ*UK* 


5^<j>  >0  <j>  ’  • 


However,  (III. 2. 36)  implies  ^hat  0  whenever 

j  6  K;  for  if  not,  then  there  would  be  an  x  e  M(F|H(A,  b),) 

such  that  X.  >  0  for  j  c  K,  contradicting  the  definition 
J 

of  K.  '  Thus  (III. 2.36) -describes  M(F1h(A,  b)),  and: 
lim  M(F|H5)  =  M(F|H(A,  b)) 

QED., 

It  would  be  easy  but  pointless  to  prove  that  if  X'^ 
solved: 


Min  /-S  log 

^<j>eJ*UK*  j 

<j>eJ*UK* 


\2>  ~  ^  V  j  e  J  U  K  , 

then  lim,  x-jf  «  x*.  where  x-^-  is  defined  in  Lemma  III. 2. 3. 
6-0+  ‘ 

There  is  a  second  way  we  can  perturb  the  problem  so 


that  it  satisfies  the  condition  that  K  be  empty.  This 
is  by  changing  the  c  values. 


Lemma  .111.2.12;  Define- new  values  for  the  Cj*s  as  follows: 

c .  if  j  e  I  U’  J 

/  c.  +  e  if  j  €  K  ' 

for  some  small  e  >  0. 

For  the  new  problem  let  the  objective  be: 

dll.2.37)  F(x)  «  S  Xj(cj  +  log  x^) 

«=  F(x)  +  es  X.  . 
jeK  J 

Then: 

(111. 2. 38)  M(FlH)  «  M(F|H),  for  each  £  >  0, 

and  if  is  the  set  of  virtual  mole  fractions  associated 
with  the  perturbed  problem,  there  exists  §(£)  e  such 
that: 

(111.2.39)  <1  V  j  (H  I  U  K. 

Note;  The  effect  of  this  perturbation  is  to  move  each  j, 
originally  in  K  to  I. 

Proof:  To  show  (111,2.38)  v/e  note  from  (III. 2. 37)  that  since 
X  >  0, 

F'(x)  =  F(x)  +  £  S  X.  >  F(x), 

jelC  J  “ 

with  equality  if  x  e  M(FlH)  (since  then  x.  =  0  if  j  e  K) . 

U6 


Thus: 


M(F|H)  “M(FlH). 

To  show  (III.  2. 39),  pick  any  5  6$*  (the  set  of  virtual 
mole  fractions  associated  with  the  unperturbed  problem) 
satisfying: 


^  for  j  €  I. 

Let  5(s)  be  defined  by: 


l^(e) 


if  j  €  I  U  J 
if  j  6  K. 


Clearly,  for  j  €  K, 


T<j>(£)  =  <  1. 

it  is  trivial  to  check  that  5(£)  e  $g.  QED. 

Again,  it  would  be  easy  but  pointless  to  prove  that 

lim,  x(t)  is  continuous  as  a  function  of  the.  ^perturbation 
t-O'*' 

e  at  e  =  0.  (The  reader  is  refered  to  Theorem  III,2i2. 
Notice  the  similarity  between  this  problem  and  that  of  prov¬ 
ing  Xg  -*  X*,  mentioned  immediately  before  Lemma  111,2.12.) 

Experience  has  led  the  author  to  believe  that  the 
solution  of  the  chemical  equilibrium  problem  is  continuous 
in  nearly  every  conceivable  perturbation,  as  one  would' 
expect  of  a  physical  system.  Unfortunately  a  proof  of  the 
general  statement  that  x(t)  x*  as  t  0  has  been  elusive. 
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